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Abstract 



We discuss the relation between the deconfining phase transition in gauge the- 
ories and the realization of the magnetic Zn symmetry. At low temperature 
the Zn symmetry is spontaneously broken while above the phase transition it 
is restored. This is intimately related to the change of behaviour of the spatial 
't Hooft loop discussed in We also point out that the realization of the 
magnetic symmetry has bearing on the behaviour of the spatial Wilson loop. 
We give a physical argument to the effect that at zero temperature the spatial 
Wilson loop must have perimeter law behaviour in the symmetric phase but 
area law behaviour in the spontaneously broken phase. At high temperature 
the argument does not hold and the restoration of magnetic Z^ is consistent 
with area law for the Wilson loop. 



1 Introduction 

This paper is devoted to further study of theoretical aspects of the deconfining 
temperature phase transition in nonabehan gauge theories. It is an immediate 
continuation of our earher work [Q. In [Q we showed that the deconfining 
phase transition in the pure Yang Mills theory is characterised by the change of 
behaviour of the 't Hooft loop operator V{C). In the "cold" phase the 't Hooft 
loop has a perimeter law behaviour < V{C) >oc exp{— aP(C)}, while in the 
"hot" phase it has an area law behaviour < V(C) >cx exp{—aS{C)}. 

In the present paper we want to sharpen somewhat this observation and further 
discuss related questions. We wish to point out that V is in fact an order 
parameter which probes the breaking of a physical symmetry of the Yang Mills 
theory. The symmetry in question is the magnetic Z]\[ symmetry discussed by 
't Hooft [Q. The deconfining phase transition is therefore characterized by the 
change in the mode of realization of a global Zjv symmetry: the symmetry is 
broken spontaneously in the "cold" phase while it is restored in the "hot" phase. 

The previous two paragraphs may sound at first like a red herring. After all 
an order parameter for the deconfining phase transition as well as a related 
Zn symmetry have been discussed for many years. The order parameter in 
question is the free energy of an external static colour source in the fundamental 
representation: the Polyakovline P = TrP exp{ig dtAo}. The symmetry 
is the transformation P exp{i^}P. We will refer to this transformation as 
the electric Zj^. There is however a great difference between the physical nature 
of P and V and the associated Zn symmetries. The operator y is a canonical 
operator in the physical Hilbert space of the Yang Mills theory. The magnetic 
Zn symmetry similarly is a transformation that acts on quantum states in the 
physical Hilbert space. On the other hand P has a very different status. It is not 
an operator in the Hilbert space and as such not a canonical order parameter. 
It appears as an auxilliary object when projecting onto gauge invariant physical 
subspace of the Hilbert space. The "electric" Z^ - the operation that transforms 
P by multiplying it by a phase - similarly is not a canonical symmetry. These 
issues were discussed in detail in |3| . There is no transformation of states in the 
physical Hilbert space that is related to this "symmetry" , although it is indeed 
a symmetry of the Euclidean path integral representing the statistical sum. 

This is not to say of course that P and electric Zjv are useless concepts. The 
standard effective action, defined by the constrained path integral 



is gauge invariant. It is instrumental in computing the vortex expectation value. 
The way the electric Z{N) symmetry is realized in Seff is also related to the 
behaviour of the order parameter of the magnetic Z{N). We will discuss this 
in detail in section 3 ptl. 
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However if one wants to describe the deconfinement phase transition in terms of 
a canonical order parameter in the same way as the Ising transition is described 
in terms of magnetisation, one should zero in on V rather than on P and should 
study the magnetic Zn symmetry rather than electric Z^. This is what we 
intend to do in this paper. 

The action of the magnetic Zjv symmetry is very different in 2+1 and 3+1 
dimensional cases. In 2+1 dimensions it acts very much like usual global sym- 
metry in a scalar theory with the order parameter being a scalar vortex field. 
In 3+1 dimensions the symmetry acts not like a standard global symetry - its 
"charge" is an integral over a two dimensional spacelike surface rather than over 
the whole of the three dimensional spaceQ As a consequence its order parame- 
ter is not a local field but rather a magnetic vortex stretching over macroscopic 
distances. 

It is therefore convenient to start the discussion with the three dimensional 
gauge theories and to present all the arguments in this case. The generalization 
of appropriate aspects of this discussion to 3+1 dimensions will be given in the 
last part of every section. 

The plan of this paper is the following. In Section 2 we recap the definition of 
the 't Hooft loop and its 2+1 dimensional analog - the magnetic vortex operator. 
We formulate the arguments for the existence of the magnetic Zjv symmetry in 
theories without fundamental matter fields. We also show by explicit construc- 
tion that the generator of this symmetry in the pure gluodynamics is none other 
than the spatial Wilson loop. 

In Section 3 wc discuss the relation between the behaviour of the 't Hooft loop 
and the realization of the magnetic Z^r in the ground state of the theory. We 
demonstrate that the mode of the realization of the symmetry changes at the 
deconfining phase transition, while spontaneously broken at low temperature 
the symmetry is restored above the phase transition. 

In Section 4 we present in a toy model a simple physical picture explaining how 
the behaviour of spatial Wilson loop discriminates at zero temperature between 
the phases with broken and unbroken magnetic Zj^. In the phase where the Zjv 
symmetry is broken, W must have an area law while in the case of unbroken 
it must have perimeter law. We explain why this argument does not generalize 
to the high temperature phase and thus why the area law behaviour of the 
Wilson loop in the hot phase is consistent with restoration of the magnetic Z]\f 
symmetry. 

Finally in Section 5 we conclude with a short discussion. 



These type of symmetries nowadays are frequently discussed in the context of "M - 
theory" M . 
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2 The magnetic symmetry and the 't Hooft 
loop operator. 



In this section we discuss the notion of the magnetic Zn symmetry and its 
order parameter - 't Hooft loop, or magnetic vortex operator. Most of the 
material contained here is not new and, perhaps with the exception of explicit 
identification of the Zjs; generator with the spatial Wilson loop, is contained in 
[H [sj H] . At the risk of being repetitive we have decided nevertheless to include 
this extended introductory part, since we feel that the concept of magnetic 
Zn symmetry is not widely appreciated in the community. The Z^r symmetry 
structure is the basis for our discussion of the deconfining phase transition in 
the following sections. 

Let us start by recalling the argument due to 't Hooft that a nonabelian SU{N) 
gauge theory with charged fields in adjoint representation posesses a global Zn 
symmetry 

We discuss the 2+1 dimensional case first. Consider a theory with several adjoint 
Higgs fields so that varying parameters in the Higgs sector the SU{N) gauge 
symmetry can be broken completely. In this phase the perturbative spectrum 
will contain the usual massive "gluons" and Higgs particles. However in addition 
to that there will be heavy stable magnetic vortices. Those are the analogs of 
Abrikosov-Nielsen-Olesen vortices in the superconductours and they must be 
stable by virtue of the following topological argument. The vortex configuration 
away from the vortex core has all the fields in the pure gauge configuration 

H°'{x) = U{x)h°', A" = iUd^U^ (2) 

Here the index a labels the scalar fields in the theory, h" are the constant vac- 
uum expectation values of these fields, and U{x) is a unitary matrix. As one 
goes around the location of the vortex in space, the matrix U winds nontrivially 
in the gauge group. This is possible, since the gauge group in the theory with- 
out fundamental fields is SU{N) /Z^^ and it has a nonvanishing first homotopy 
group IIi{SU{N)/Zn) = Z^- Practically it means that when going around the 
vortex location full circle, U does not return to the same SU{N) group element 
Uq, but rather ends up at exp{i^}Uo. Adjoint fields do not feel this type of 
discontinuity in U and therefore the energy of such a configuration is finite. 
Since such a configuration can not be smoothly deformed into a trivial one, a 
single vortex is stable. Processes involving annihilation of N such vortices into 
vacuum are allowed since N-vortex configurations are topologically trivial. One 
can of course find explicit vortex solutions once the Higgs potential is specified. 
As any other semiclassical solution in the weak coupling limit the energy of such 
a vortex is inversely proportional to the gauge coupling constant and therefore 
very large. One is therefore in a situation where the spectrum of the theory 
contains a stable particle even though its mass is much higher than masses of 
many other particles (gauge and Higgs bosons) and the phase space for its decay 
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into these particles is enormous. The only possible reason for the existence of 
such a heavy stable particle is that it must carry a conserved quantum number. 
The theory therefore must possess a global symmetry which is unbroken in the 
completely higgsed phase. The symmetry group must be Zj^ since the number 
of vortices is only conserved modulo N. 

Now imagine changing smoothly the parameters in the Higgs sector so that the 
expectation values of the Higgs fields become smaller and smaller, and finally the 
theory undergoes a phase transition into the confining phase. One can further 
change the parameters so that the adjoint scalars become heavy and eventually 
decouple completely from the glue. This limiting process does not change the 
topology of the gauge group and therefore does not change the symmetry content 
of the theory. We conclude that the pure Yang-Mills theory also posesses a Zn 
symmetry. Of course since the confining phase is separated from the completely 
Higgsed phase by a phase transition one may expect that the Zjv symmetry in 
the confining phase is represented differently. In fact the original paper of 't 
Hooft as well as subsequent workQ convincingly argued that in the confining 
phase the Z^ symmetry is spontaneously broken and this breaking is related to 
the confinement phenomenon. 

The physical considerations given above can be put on firmer formal basis. In 
particular one can construct explicitly the generator of the Zj^ as well as the 
order parameter associated with it - the operator that creates the magnetic 
vortex Iql. We will now describe this construction. 



2.1 The Abelian case 

Consider first an Abelian gauge theory. In this case the homotopy group is 
Z and therefore we expect the U{1) rather than Zpf magnetic symmetry. It 
is in fact absolutely straightforward to identify the relevant charge. It is none 
other than the magnetic flux through the equal time plane, with the associated 
conserved current being the dual of the electromagnetic field strength 

$ = y" d^xB{x), d^Ff, = (3) 

The current conservation is insured by the Bianchi identity. A group element 
of the C/(l) magnetic symmetry group is exp{ia<I>} for any value of a. A local 
order parameter - a local field V{x) which carries the magnetic charge - is also 
readily constructed. It has a form of the singular gauge transformation operator 
with the singularity at the point x 



V{x) — exp - / (fy 
9 J 



E,{y) + Q{x~y)My) 



(4) 



where Q{x — y) is the polar angle function and Jq is the electric charge density 
of whatever matter fields are present in the theory. The cut discontinuity in 
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the function Q is not physical and can be chosen parallel to the horizontal axis. 
Using the Gauss' law constraint this can be cast in a different form, which we 
will find more convenient for our discussion 

V{x) = exp / dy' e^jE,{y) (5) 

9 Jc 

where the integration goes along the cut of the function which starts at 
the point x and goes to spatial infinity. The operator does not depend on 
where precisely one chooses the cut to lie. To see this, note that changing the 
position of the cut C to C" adds to the phase ^ dPxdiE^ where S is the 
area bounded by C — C". In the theory we consider only charged particles with 
charges multiples of g are present. Therefore the charge within any closed area 
is a multiple integer of the gauge coupling cPxdiE^ — gn and the extra phase 
factor is always unity. 

The meaning of the operator V is very simple. From the commutation relation 

Vix)B{y)V\x)^Biy) + jS'ix-y) (6) 

it is obvious that V creates a pointlike magnetic vortex of flux 27r/g. Despite its 
nonlocal appearance the operator V can be proven to be a local Lorentz scalar 
field [p^. The locality is the consequence of the fact that V{x) commutes with 
any local gauge invariant operator in the theory 0{y) except when x = y. This is 
due to the coefficient 27r/(/ in the exponential which ensures that the Aharonov- 
Bohm phase of the vortex created by V and any dynamical charged particle 
present in the theory vanishes. Eqs.(||j|) formalize the physical arguments of 't 
Hooft in the abclian case. 

2.2 The non-Abelian case at weak coupling. 

Let us now move onto the analogous construction for nonabelian theories. Ul- 
timately we are interested in the pure Yang - Mills theory. It is however il- 
luminating to start with the theory with an adjoint Higgs field and take the 
decoupling limit explicitly later. For simplicity we discuss the SU (2) gauge the- 
ory. Consider the Georgi-Glashow model - SU (2) gauge theory with an adjoint 
Higgs field. 

where 

VfH^ = d^H'' ~ gf^^AlH" (8) 

At large and positive jl^ the model is weakly coupled. The SU{2) gauge sym- 
metry is broken down to U{1) and the Higgs mechanism takes place. Two gauge 
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bosons, W^, acquire a mass, while the third one, the "photon", remains mass- 
less to all orders in perturbation theory. The theory in this region of parameter 
space resembles very much electrodynamics with vector charged fields. The 
Abelian construction can therefore be repeated. The SU{2) gauge invariant 
analog of the conserved dual field strength is 

= ^[e^-AF^A"" - -e^'^^e^^'^naiV^nfiVxriy] (9) 
^ 9 

where = -j^ is the unit vector in the direction of the Higgs field. Classically 
this current satisfies the conservation equation 

d^F^ = (10) 

The easiest way to see this is to choose a unitary gauge of the form H°-{x) = 
H{x)S'^^. In this gauge F is equal to the abehan part of the dual field strength 
in the third direction in colour space. Its conservation then follows by the 
Bianchi identity. Thus classically the theory has a conserved U{1) magnetic 
charge $ = / cPxFq just like QED. However the unitary gauge can not be 
imposed at the points where H vanishes, which necessarilly happens in the 
core of an 't Hooft-Polyakov monopole. It is well known of course |§] that the 
monopoles are the most important nonperturbative configurations in this model. 
Their presence leads to a nonvanishing small mass for the photon as well as to 
confinement of the charged gauge bosons with a tiny nonperturbative string 
tension. As far as the monopole effects on the magnetic fiux, their presence 



leads to a quantum anomaly in the conservation equation (10). As a result only 
the discrete Z2 subgroup of the transformation group generated by $ remains 
unbroken in the quantum theory. The detailed discussion of this anomaly, the 
residual Z2 symmetry and their relation to monopoles is given in |6|. 

The order parameter for the magnetic Z2 symmetry is constructed analogously 
to QED as a singular gauge transformation generated by the gauge invariant 
electric charge operator 

= ef^'^^d.iF^n''), Q = J d^xJo{x) (11) 



Explicitly 



V(x) = exp - / cPy 
9 J 



+ e(a; - y)Uy) 

[x-yY 



= exp / dfe^.n^'E^iy) (12) 

9 Jc 

One can think of it as a singular SU{2) gauge transformation with the field 
dependent gauge function 

X-{y)^-e{x-y)n-{y) (13) 
9 
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This field dependence of the gauge function ensures the gauge invariance of 
the operator V. Just like in QED it can be shown[H, that the operator 
y is a local scalar field. Again like in QED, the vortex operator y is a local 
eigenoperator of the abelian magnetic field B{x) — Fq. 

[V{x),B{y)]^~—V{x)6^ix-y) (14) 
9 

That is to say, when acting on a state it creates a pointlike magnetic vortex 
which carries a quantized unit of magnetic flux. The Z2 magnetic symmetry 
transformation is generated by the operator 

U = exp{i|$} (15) 

and acts on the vortex field ^ as a phase rotation by n 

e'^*F(a;)e-^** = -V{x) (16) 



An operator closely related to U and which will be of interest to us in the 
following, is the generator of the magnetic Z2 transformation only inside some 
closed contour C 

U{C) = exp{i^ / d?xB{x)) (17) 
2 Js 

where the integration is over the area S bounded by C. The analog of the 
commutator eq.([T^) for this operator is 

UcV{x)ul = - V{x) , xeS 

V{x) , xiS (18) 

Taking the contour C to run at infinity Uc becomes the generator of Z2. 

We now have the explicit realization of the magnetic Z2 symmetry in the Gcorgi- 
Glashow model. 



2.3 The pure gauge theory. 

Our next step is to move on to the pure Yang Mills theory. This is achieved 
by smoothly varying the /i^ coefficient in the Lagrangian so that the coeficient 
of the mass term of the Higgs field becomes positive and eventually arbitrarily 
large. It is well known that in this model the weakly coupled Higgs regime and 
strongly coupled confining regime are not separated by a phase transition|Q. 
The pure Yang Mills limit in this model is therefore smooth. 

In the pure Yang Mills limit the expressions eq. (^,|l^,^7|) have to be taken with 
care. When the mass of the Higgs field is very large, the configurations that 
dominate the path integral are those with very small value of the modulus of 
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the Higgs field \H\ cx 1/M. The modulus of the Higgs field in turn controls the 
fluctuations of the unit vector n", since the kinetic term for n in the Lagrangian 
is |i?p(D^n)^. Thus as the mass of the Higgs field increases the fluctuations of 
n grow in both, amplitude and frequency and the magnetic field operator B as 
defined in eq.(^) fluctuates wildly. This situation is of course not unusual. It 
happens whenever one wants to consider in the effective low energy theory an 
operator which explicitly depends on fast, high energy variables. The standard 
way to deal with it is to integrate over the fast variables. There could be two 
possible outcomes of this procedure. Either the operator in question becomes 
trivial (if it depends strongly on the fast variables), or its reduced version is 
well defined and regular on the low energy Hilbert space. The "magnetic field" 
operator B in eq.(^ is obviously of the first type. Since in the pure Yang 
Mills limit all the orientations of n° are equally probable, integrating over the 
Higgs field at fixed will lead to vanishing of B. However what interests us 
is not so much the magnetic field but rather the generator of the magnetic Z2 
transformation Uc of eg. (|T7|) . In the pure Yang-Mills limit we are thus lead to 
consider the operator 

Uc = liniH^a J Dn^ exp | - |ff|2(^n,)2| (19) 

exp {i| £ d^x{e,,Ftjn- - ^e'^ e'''^na{V,n)\V,ny)} 

The weight for the integration over n is the kinetic term for the isovector Ha- 
As was noted before the action does not depend on n°- in the YM limit. This 
term however regulates the integral and we keep it for this reason. This op- 
erator may look somewhat unfamiliar at first sight. However in a remarkable 
paper jl^ Diakonov and Petrov showed that eq.(20) is equal to the trace of the 
fundamental Wilson loop along the contour C^. 

Uc^Wc = TtV exp |ig y" dl'A'l (20) 



We conclude, that in the pure Yang-Mills theory the generator of the magnetic 
Z2 symmetry is the fundamental spatial Wilson loop along the boundary of the 
spatial plane. 

There is a slight subtlety here that may be worth mentioning. The generator 
of a unitary transformation should be a unitary operator. The trace of the fun- 
damental Wilson loop on the other hand is not unitary. One should therefore 
strictly speaking consider instead a unitarizcd Wilson loop W = -^M—. How- 
ever the factor between the two operators vWW^ is an operator that is only 
sensitive to behaviour of the fields at infinity. It commutes with all physical 

^We note that Dyakonov and Petrov had to introduce a regulator to defin e th e path integral 
over n. The regulator they required was precisely of the same form as in eq. (bol). It is pleasing 
to see that this regulator appears naturally in our approach as the remnant ofthe kinetic term 
of the Higgs field. 
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local operators 0{x) unless x oo. In this it is very different from the Wilson 
loop itself, which has a nontrivial commutator with vortex operators V{x) at all 
values of x. Since the correlators of all gauge invariant local fields in the pure 
Yang Mills theory are massive and therefore short range, the operator vWW^ 
must be a constant operator on all finite energy states. The difference between 
W and W is therefore a trivial constant factor and we will not bother with it in 
the following. Perhaps of more concern is the difference between Wc and Wq 
when the contour C is not at infinity. However here again the factor between 



the two operators yWcW^. is only sensitive to physical degrees of freedom on 
the contour C and not inside it. Due to its presence the vacuum averages of 
Wc and Wc may differ at most by a factor which has a perimeter behaviour 
< Wc >= exp{mP(C)} < Wc > where P(C) is a perimeter of C. The question 
we will be interested in is whether < Wc > has a perimeter or area behaviour. 
As far as the answer to this question is concerned Wc and Wc are completely 
equivalent, and we will not make distinction between them. In the rest of this 
paper we will therefore refer to W as the generator of Z2 keeping this little 
caveat in mind. 



Next we consider the vortex operator eq.(^J|). Again we have to integrate it 
over the orientations of the unit vector n". This integration in fact is equivalent 



to averaging over the gauge group. Following |11 one can write Ua in terms of 
the SU(2) gauge transformation matrix Q. 

n = ^Tinrn^Ts (21) 
The vortex operator in the pure gluodynamics limit then becomes 

V{x)= [ DO exp — / dyie^jTiflE^n^Ts (22) 
J 9 Jc 

This form makes it explicit that V{x) is defined as the gauge singlet part of the 
following, apparently non gauge invariant operator 

V{x)=exp— dfe,,Ef{y) (23) 
9 Jc 

The integration over obviously projects out the gauge singlet part of V. In the 
present case however this projection is redundant. This is because even though 
V itself is not gauge invariant, when acting on a physical state it transforms it 
into another physical state|^. By physical states we mean the states which satisfy 
the Gauss' constraint in the pure Yang-Mills theory. This property of V was 
noticed by 't Hooft ||]. To show this let us consider V{x) as defined in eq.(p3|) 
and its gauge transform = ^l^VH, where ft is an arbitrary nonsingular gauge 
transformation operator. The wave functional of any physical state depends 

■^This is not a trivial statement, since a generic nongauge invariant operator has nonvan- 
ishing matrix elements between the physical and an unphysical sectors. 
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only on gauge invariant characteristics of the vector potential, i.e. only on the 
values of Wilson loops over all possible contours. 

MA^] = mW{C)}] (24) 

Acting on this state by the operators V and Vh respectively we obtain 

V\^> = ^v[Ai] = 'f[{VW{C)V^}] 
Vn\-^> = ¥^[A,]^mVnW{C)V^,}] (25) 

It is however easy to see that the action of V{x) and Vn{x) on the Wilson 
loop is identical - they both multiply it by the centergroup phase (which stays 
unaffected by fl) if x is inside C and do nothing otherwise. Therefore we see 
that 

]/|*>=yo|*> (26) 
for any physical state 5*. Thus we have 

>= f^yf^^lf >= yi* > (27) 

where the first equality follows from the fact that a physical state is invariant 
under action of any gauge transformation H. and the second equality follows 
from eq.(p6|). But this equation is nothing but the statement that the state 
yjvl/ > is physical, i.e. invariant under any nonsingular gauge transformation. 

We have therefore proved that when acting on a physical state the vortex op- 
erator creates another physical state. For an operator of this type the gauge 
invariant projection only affects its matrix elements between unphysical states. 
Since we are only interested in calculating correlators of V between physical 
states, the gauge projection is redundant and we can freely use V rather than 
V to represent the vortex operator. 

It is instructive to note that this propery is not shared by the Wilson loop. One 
can in fact represent the Wilson loop as a singlet gauge projection of a simple 
Abelian loop operator. The second exponential in eq.(pO[) can be written as 

exp{z|y" dPAln'^ - ^ J d^xe,je'''''nad,nbdjn,} (28) 

Using eq.(|2l]) we can rewrite the integral in eq.(|2^) -omitting the regulating 
kinetic piece- as: 

Wc^ I -OrJexpl^l I TrTz{nA'-n'^ + ind'n'^)dU^ (29) 



The Wilson loop is therefore the gauge singlet part of the Abelian loop 

= expi| / TrAVadP (30) 
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The matrix elements of Wc and Uq on physical subspace therefore are the same. 
However Uq as opposed to V does have nonvanishing nondiagonal matrix ele- 
ments, that is matrix elements between the physical and the unphysical sectors. 
It therefore can not be used instead of Wc in gauge theory calculations. For 
example non gauge invariant states will contribute as intermediate states in the 
calculation of quantities like the correlation function < Uq^Uq^ >, while their 
contribution vanishes in similar correlators which involve the Wilson loop. 

The generalization of the preceding discussion to SU{N) gauge theories is 
straightforward. Once again one can start with the Georgi-Glashow like model, 
where the SU{N) is higgsed to ?7(1)*^^~^^0. The construction of the vortex 
operator and the generator of Zj^ in this case is very similar and the details are 
given m §. Taking the mass of the Higgs field to infinity again projects the 
generator onto the trace of the fundamental Wilson loop. The vortex operator 
can be taken as 

47ri /■ 

V{x) = exp{— d2/^e,,Tr(YEi(y)) (31) 

where the hypercharge generator Y is defined as 

r = diag(l,l,...,-(7V-l)) (32) 

and the electric field is taken in the matrix notation Ei — X'^Ef with A" - the 
SU (N) generator matrices in the fundamental representation. 



2.4 Generalization to 3+1 dimensions 

To conclude this section we discuss how the magnetic symmetry structure gen- 
eralizes to four dimensions. The conserved Zn generator in the Georgi-Glashow 
model is defined through 

Us = exp{^|^d2^^^B>" - ^e'^''e'''^na{Vjn)\Vkny)} (33) 

Although the definition of U contains explicitly the surface S through which 
the abelian magnetic fiux is integrated, the operator in fact does not depend 
on S but is specified completely by its boundary. This is because changing S 
changes the phase of U by the magnetic flux through the closed surface. The 
only dynamical objects that carry magnetic flux in the theory are 't Hooft- 
Polyakov monopoles. Since their flux is quantized in units of 4:TT/g the change 
in the phase is always a multiple integer of 27r. In the pure Yang-Mills limit 
the operator Us again reduces to the trace of the fundamental Wilson loop 

*In SU (N) theories with N > 2 there in principle can be phases separated from each other 
due to spontaneous breaking of some global symmetries. For instance the SU (3) gauge th eory 
with adjoint matter has a phase with spontaneously broken charge conjugation invariance |17[ . 
Still even in this phase the confining properties are the same as in the strongly coupled pure 
Yang-Mills theory, with the Wilson loop having an area law. 
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along the boundary of S. Taking the contour to infinity defines the generator 
of magnetic Zff. As we have aheady noted, this charge is a little unusual in 
that it is defined as a surface integral rather than a volume integral. As a result 
the order parameter for this symmetry transformation is not a local but rather 
a stringy field. This is of course just a restatement of the fact that magnetic 
vortices in 3+1 dimensions are stringlike objects. The operator that creates a 
vortex can still be defined in a way similar to 2+1 dimensions. Skipping the 
intermediate steps which we went through in the previous discussion we give the 
final result for the pure Yang Mills SU (N) gauge theory. The magnetic vortex 
along the curve C is created by the following operator of the "singular gauge 
transformation" 

V{C) = exp{-^ J d'xTT(D'u;cYW} = exp{^ ^ d2s'Tr(YE')} (34) 

with ujc{x)j the singular gauge function which is equal to the solid angle sub- 
tended by C as seen from the point x. The function lo is continuos everywhere, 
except on a surface S bounded by C, where it jumps by 47r. Other than the 
fact that S is bounded by C, its location is arbitrary. The vortex loop and the 
spatial Wilson loop satisfy the 't Hooft algebra 

V^{C)WiC')V{C) = e^"('^'^')Ty(C") (35) 

where n{C, C) is the linking number of the curves C and C". One can consider 
closed contours C or infinite contours that run through the whole system. For 
an infinite contour C and the Wilson loop along the spatial boundary of the 
system the linking number is always unity. The V{C) for an infinite loop is 
therefore an eigenoperator of the Zn magnetic symmetry and is the analog of 
the vortex operator V{x) in 2+1 dimensions. Any closed vortex loop of fixed 
size commutes with the Wilson loop if the contour C is very large. Such a 
closed loop is thus an analog of the vortex-antivortex correlator V{x)V'^ {y), 
which also commutes with the global symmetry generator, but has a nontrivial 
commutator with Uc if C encloses only one of the points x or y. 

To summarize this section, we have shown that pure Yang Mills theory in 2+1 
and 3+1 dimensions has a global Zn magnetic symmetry. The generator of the 
symmetry group in both cases is the trace of the fundamental Wilson loop along 
the spatial boundary of the system. The order parameter for this symmetry in 
2+1 dimensions is a local scalar field V{x), while in 3+1 dimensions a stringlike 
field V{C). In both cases the field V is gauge invariant on physical states and is 
a bona fide canonical order parameter which distinguish in gauge invariant way 
the phases of the theory. In the next section we discuss the realization of the 
magnetic symmetry in the confined and the deconfined phases. 

^The derivative term d^ui in this expression should bo understood to contain only the 
smooth part of the derivative and to exclude the contribution due to the discontinuity of uj 
on the surface S. 
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3 Hot and cold realizations of the magnetic Zjy. 



As with any global symmetry, it is important to understand what is the mode 
of realization of magnetic in the ground state of the theory. This mode of 
realization depends of course on the parameters of the theory as well as on the 
temperature. The situation at zero temperature is well understood. 

3.1 2+1 dimensions. 

Again we start with three dimensions. There is a very general argument due to 
't Hooftj2||^ stating that if the theory does not have zero mass excitations the 
area law of the Wilson loop implies the nonvanishing expectation value of the 
vortex operator V{x). Conversely if the Wilson loop has a perimeter law the 
expectation value of V{x) must vanish and the correlation function V {x)V^ {y) 
must have an exponential falloff with |a; — It follows that in the pure Yang 
Mills theory the vacuum expectation value of the vortex operator does not vanish 
and therefore the Z]\[ magnetic symmetry is spontaneously broken. The same 
is true in the partially broken Higgs phase of the Georgi-Glashow model. As 
mentioned in the last section the confining and the Higgs regimes in this model 
are analytically connected and therefore the realization of all global symmetries 
in the two regimes is the same. 

In fact in the weakly coupled Higgs phase this can be verified by the direct calcu- 
lation of the expectation value of y Q . This calculation maps very simply into 
the classic monopole plasma calculation of Polyakov and was discussed in detail 
in 1^. One can also explicitly construct the low energy effective Lagrangian in 
terms of the field V which realizes the spontaneously broken Zj^ symmetry and 
describes the low energy spectrum of the Georgi - Glashow vacuum. 

C ^ d^V*d''V - X{V*V - - CiV^ + V*^) (36) 

A similar effective Lagrangian with some quantitative differences was argued to 
be valid also for the pure Yang-Mills theory in p^ . 

The application of the 't Hooft argument at finite temperature is somewhat 
less straightforward. Since at finite temperature the Lorentz invariance is bro- 
ken, the temporal and spatial Wilson loops do not necessarily have the same 
behaviour and one has to be more careful. The original argument relates the 
behaviour of the vortex operator and the temporal Wilson loop. At finite tem- 
perature in the Euclidean formalism the extent of the system in the temporal 
direction is finite. As a result it is not possible to distinguish between the area 
and perimeter law for "asymptotically" large temporal loops. Instead the role 
of the temporal Wilson loop is taken over by the Polyakov line - the loop that 

^The original argument as stated in is formulated for 3+1 dimensional theories, however 
its generalization to 2+1 dimensions requires only linguistic changes. 
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winds around the total volume of the system in the temporal direction. Thus 
one expects that in the deconfining phase where the Polyakov line has a nonva- 
nishing vacuum average, the vortex operator should have vanishing expectation 
value. Indeed this can be easily confirmed by the explicit calculation of the 
VEV of the vortex operator using the method of . In the calculation was 
performed in 3+1 dimensions, but adapting it to 2+1 dimensional case is trivial. 
We give below a brief outline. 

Consider the equal time vortex-antivortex correlation function. At finite tem- 
perature is it given by the following expression 

< V{x)VHy) >= Tre-^(^'+^')e^'^/>-''''^^ (37) 

The line integral can be taken along the straight line L connecting the points 
X and y. For definiteness we take x and y to be separated in the direction of 
the first axis. Introducing the imaginary time axis and the Lagrange multiplier 
field Aq in the standard way this expression can be transformed to 

<V{x)VHy) >= J DADAoexp{-^ dt j Sxi^d^AI - {D^A^T - 5(t)aff+{B'^f} 

(38) 

where the "external field " a' is given by 

<(x)=<5''3^,2— '5(x-X) (39) 
9 

The delta function in time in front of the external field a* in eq. (^) appears for 
the following reason. As we saw in section 2.3 the product of the vortex and 
an antivortex operator is gauge invariant. This is because it induces a singular 
gauge transformation which is continuous up to the center element. However if 

we split it up in imaginary time into infinitesimal bits Vat = e*''*"^ ' '^'^'^^ 
then any single such bit separately is not gauge invariant, since the transforma- 
tion it induces is genuinely discontinuous across the line connecting the points 
X and y. The operators Vdt therefore do not commute with the projection op- 
erator on physical states. To obtain the path integral representation for the 
expectation value eq.(^7|) we should introduce the projection operator only at 
the last point in imaginary time and not at the intermediate points. In the 
path integral language this corresponds to the gauge fixing Aq = Q everywhere 
except at one time, say t = 0. In this gauge it is straightforward to see that 
the Gaussian integration over the electric field leads to the usual path integral 
with diAt) shifted by a^. This can then be rotated to an arbitrary gauge with 
the result eq.(||)|. 

^In this derivation we dropped commutator terms between the Hamiltonian and the ex- 
ponent in the vortex-antivortex operator. These commutator terms only exist at t = and 
therefore drop out in the continuum hmit (i.e are 0{dt)). In the lattice realization they are 
indeed present and complete the expression eq.rtSSl) to the twisted plaquette representation [W 



15 



To evaluate the path integral eq.(|3^) we follow the standard procedure and inte- 
grate out all modes except for the Polyakov loop in a saddle point approximation 
This leads to the effective action Seff{q,ai), where \/2TrP — cosq 

< V{x)V\y) J Dqex_p-Seff{q,a,) (40) 

To one loop order the effective action is given by 

5e// = J <fx{^{d,q+^-a,f + U{q)) (41) 

The matrices t°- in are the generators of SU{2) in the fundamental repre- 
sentation and are normalized according to trr'^r'' — ^6°'''. 



The one loop effective potential U is related to a Bernoulli polynomial and 
can be read off the expressions in iQ. The only property of U which is 
important to us is that it has two degenerate minima at g = 0, tt. 

To calculate the correlator we have to find the configuration of q which mini- 
mizes the action eg. (pl|) . Qualitatively the form of the solution is clear. The 
considerations identical to those in ||l| tell us that it must be the "broken" elec- 
tric Z2 domain wall: half a wall ( q -^x2^co 0, q{x2 = 0) = -1) above the line L 
and half a wall ( q{x2 = 0) = — -j, q -^x^^-oo 0,) below the line L separated by 
a discontinuity Sq = n. The action of such a configuration is Seff = a\x — y\ 
where a is the domain wall tension" . The vortex correlator is thus given 

by 

< V{x)V''{y) >= exp{-a|x - y|} (42) 

As |x— y I become large the correlation function decreases exponentially, and thus 
the expectation value of the vortex operator vanishes. For the SU{N) group 
this calculation trivially generalizes and gives the same result. The exponential 
decay is also obtained for the correlator of V"^ with any power m < N. 

Recall that the vortex operator is the order parameter for the magnetic Zn sym- 
metry. Moreover the powers of V exhaust all possible local order parameter^. 
Their vanishing is therefore an unambiguous indication that the magnetic Zn 
is restored in the high temperature deconfined phase. 

In hindsight this is not very surprising. Indeed, we are dealing with physical 
discrete symmetry which is spontaneously broken at zero temperature. When 
the system is heated it is unavoidable that entropy effects take over and at 
some sufficiently high temperature the symmetry must be restored. A good 



qualitative guide here is the effective Lagrangian eq.(36). It describes a simple 
Z2 invariant scalar theory. There is very little doubt that a system described by 
this Lagrangian indeed undergoes a symmetry restoring phase transition at some 

*The latter statement is correct modulo multiplication of V" by local gauge invariant and 
Zpf invariant operators. These possible factors do not change the fact of the exponential decay 
of the corelators and are therefore unimportant for our discussion. 
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Tc- Moreover the effective Lagrangian approach also suggests that this phase 
transition has deconfining character. As shown in |lO[ ^ the charged states 
in the effective theory eq. (|3^) are represented by sohtonic configurations of the 
vortex field V with unit winding number. The energy of any such state is linearly 
divergent in the infrared. The reason is that due to finite degeneracy of vacuum 
states, the minimum energy configuration looks like a quasi onedimensional strip 
across which the phase of V winds. The energy density inside this "electric 
flux tube" is proportional to the vacuum expectation value of V. When the 
VEV vanishes, so does the string tension. Stated in other words, when V 
vanishes, the phase fluctuations are large and the winding number is not a 
sharp observable. An external charge is thus screened easily by regions of space 
around it with vanishing V . The phase with < V >= is therefore not 
confining. In the theory with several Higgs fields this phase exists even at zero 
temperature and corresponds to a completely Higgsed phase - where the gauge 
group is broken completely. In such a Higgs phase indeed the colour is screened 
rather than confined. In the pure Yang-Mills theory this phase is absent at zero 
temperature, but is realized as the deconfined phase at T > Tc- We thus see that 
the behaviour of the vortex operator at high temperature does indeed match 
the simple intuition coming from a invariant effective Lagrangian very well. 

3.2 Extension to 3+1 dimensions 

The 't Hooft argument now states that the vanishing vacuum average of the 
Polyakov line is incompatible with the area law behaviour of the spatial 't Hooft 
loop and vice versa. This means that in the confining phase the 't Hooft loop 
has perimeter law. In the high temperature deconfined phase the behaviour of 
the spatial 't Hooft loop must become area since the average of the Polyakov 
line is finite. Again this is confirmed by explicit calculation in Q. 

A more subtle question is how the behaviour of the 't Hooft loop relates to the 
realization of the magnetic Zn symmetry. The Zn symmetry does not have 
an order parameter which is a local field defined at a point. The only order 
parameters in the strict sense (an eigenoperator with a nonvanishing eigenvalue) 
is a 't Hooft line V{C) which runs through the whole system p^ . 

In a system which is finite in the direction of the loop, but is infinite in the 
perpendicular directions everything is clear cut. In this case there are two 
possibilities: 

a) < >7^ and the magnetic Zn broken, or 

b) . < V >— and the nagnetic Z^ restored. 

In the system infinite in all directions C is necessarily an infinite line, and the 
expectation value < V{C) > clearly vanishes irrespective of whether Z^ is 
broken or not. The 't Hooft loop along a closed contour on the other hand 
is never zero, since it is globally invariant under the Zj\j transformation. It is 
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therfore impossible to find an operator whose VEV distinguishes between the 
two possible realizations of the magnetic symmetry by vanishing in one phase 
and not vanishing in the other. Nevertheless the behaviour of the closed loop 
does indeed reflect the realization mode of the symmetry, since it is qualitatively 
different in the two possible phases. Namely vacuum expectation value of a large 
closed 't Hooft loop (by large, as usual we mean that the linear dimensions of 
the loop are much larger than the correlation length in the theory) has an area 
law decay if the magnetic symmetry is spontaneously broken, and perimeter law 
decay if the vacuum state is invariant. 

To understand the physics of this behaviour it is useful to think of the 't Hooft 
line as built of "local" operators - httle "magnetic dipoles". Consider eq.(|^ 
with the contour C running along the x- axis and the surface S chosen as the 
(a;, y) plane. Let us mentally divide the line into (short) segments of length 
2 A centered at Xi. Each one of these segments is a little magnetic dipole and 
the 't Hooft loop is a product of the operators that create these dipoles. The 
definition of these little dipole operators is somewhat ambiguous but since we 
only intend to use them here for the purpose of an intuitive argument any 
reasonable definition will do. It is convenient to define a single dipole operator 
in the following way 

Da{x) = exp{i J d^'yia+ix + A ~ y) + ai{x - A ~ y)]Tr(YE'(y))} (43) 

where a^{x — y) is the c-number vector potential of the abclian magnetic 
monopole (antimonopole) of strength Air/gN . The monopole field correpsonding 
to Gi contains both, the smooth Xi jx^ part as well as the Dirac string contribu- 
tion. The Dirac string of the monopole - antimonopole pair in eq. (^) is chosen 
so that is connects the points x — A and a; + A along the straight line. The 
dipole operators obviously have the property 

DA(a;)L'A(a; + 2A) =i?2A(a; + A) (44) 

This is because in the product the smooth field contribution of the monopole in 
Di^{x) cancels the antimonopole contribution in Da{x + 2 A), while the Dirac 
string now stretches between the points {x — A) and {x + 3A) . When multiplied 
over the closed contour, the smooth fields cancel out completely, while the sur- 
viving Dirac string is precisely the magnetic vortex created by a closed 't Hooft 
loop operator. The 't Hooft loop can therefore be written as 

V (C) =U,^D Aix^) (45) 

The dipole operator D{xi) is an eigenoperator of the magnetic flux defined 
on a surface that crosses the segment [xi — A,Xi + A]. Suppose the magnetic 
symmetry is broken. Then we expect the dipole operator to have a nonvanishing 
expectation valuej^ < D >= d{A). If there are no massless excitations in the 

^The magnetic dipole operators defined above are strictly speaking not local, since they 



18 



theory, the operators D{xi) and D{xj) should be decorrelated if the distance 
Xi — Xj is greater than the correlation length I. Due to eq.(|45|), the expectation 
value of the 't Hooft loop should therefore roughly behave as 



<y(C)>=d(0^/'=exp{-ln(^— jy} (46) 

where L is the perimeter of the loop. In the system of finite length Lx, the 
vacuum expectation value of the vortex line which winds around the system in 
x-direction is therefore finite as in eq.(p6|) with L ^ L^- 



On the other hand in the unbroken phase the VEV of the dipole operator de- 
pends on the size of the system in the perpendicular plane Ly. For large Ly it 
must vanish exponentially as d = exp{— aLj,}. So the expectation value of V 
behaves at finite Ly in the unbroken phase as: 

< ViC) cxp{-aL,^Lj (47) 

and vanishes as — > oo. Thus in a system which is finite in x direction, but 
infinite in y direction, the 't Hooft line in the x direction has a finite VEV in 
the broken phase and vanishing VEV in the unbroken phase. 

In the limit of the infinite system size oo the VEV obviously vanishes in 

both phases. This is of course due to the fact that y is a product of infinite 
number of dipole operators, and this product vanishes even if individual dipole 
operators have finite VEV^. However one can avoid any reference to finite size 
system and infinite vortex lines by considering closed 't Hooft loops. 

For a closed loop with long sides along x axis at y = and y = R the above argu- 
ment leads to the conclusion that in the broken phase V must have a perimeter 
law, eq.(|4^). In the unbroken phase the correlation between the dipoles at y = 
and dipoles at y — R should decay exponentially < D{0)D{R) >(x exp{— aj} 
and thus 

LR^ , 
-jT} = exp{-a- 

Thus the perimeter behaviour of < V{C) > indicates a vacuum state which 
breaks spontaneously the magnetic while the area behaviour means that 
the magnetic Zn is unbroken. 

The results of then mean that in 3-1-1 dimensions as well as in 2-1-1 dimension 
the magnetic symmetry is restored above the deconfining phase transition, in 
the sense of eq.( |47|). 

In the next section we discuss what is the implication of this conclusion on the 
behaviour of the spatial Wilson loop. 

carry the long range magnetic field of a dipole. However, the dipole field falls off with distance 
very fast. Therefore even though this fall off is not exponential the slight nonlocality of D 
should not affect the following qualitative discussion. 

^''The VEV of the dipole D must be smaller than one since D is defined as a unitary 
operator. 



< V{C) exp{-a— } = exp{-a-} (48) 
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4 Spatial Wilson loop at low and high temper- 
ature. 



As we have shown in Section 2 the spatial Wilson loop is the generator of the 
magnetic Zn symmetry. We expect therefore that the mode of realization of 
the magnetic Zn is strongly linked to the behaviour of W. The argument is 
simplest to state for a toy model which exemplifies the basic physics in a very 
simple setting. 

Rather than talk about nonabelian gauge theory, consider a scalar theory of a 
complex field ^ with global Zn theory in 2+1 dimensions. 

C = d^(l>d^r + A(<^> - ii'^f + C(<A^ + (</)*)^) (49) 

The generator of the Ziq symmetry is given by 

U = exp {i^ j dPxjoix)] = exp y d['x{'Kcj> - it*4>*)] (50) 

where tt = do(p* is the momentum conjugate to the field (j). Obviously with the 
canonical commutation relations between tt and (j) one has 

U(l){x)U'^ =e'^(l){x) (51) 

We will be interested in the behaviour of the operator which generates the Zjv 
transformation only inside some region S of the two dimensional plane. 

U{S) = exp[^ f d^xiTT<P-7T*<P*)} (52) 
s 

U{S)^{x)U^S) = e'^(p{x) xgS 
= dix), x^S 

We will refer to this operator as the U-loop. Throughout this discussion we 
assume that there are no massless excitations in the spectrum of the theory and 
that the linear dimensions of the area S are much larger than the correlation 
length. 

The statement we are aiming at is that at zero temperature in the phase with 
broken Zi^ the U-loop has an area law behaviour while in the phase with un- 
broken Zjv this changes into the perimeter law behaviour. 

4.1 U-loop in the broken phase 

Consider the broken phase first. Wc arc interested in the vacuum expectation 
value of U{S). This is nothing but the overlap of the vacuum state < 0| and the 
state which is obtained by acting with U {S) on the vacuum state >= C/|0 >. 
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If the symmetry is broken, the field 4> in the vacuum state is pointing in some 
fixed direction in the internal space. In the state jS* > on the other hand its 
direction in the internal space is different - rotated by 27r/iV - at points inside 
the area 5'. In the local theory with finite correlation length the overlap between 
the two states approximatelly factorises into the product of the overlaps taken 
over the region of space of linear dimension of order of the correlation length I 

<0|5>-n, <0,|5, > (53) 

where the label x is the coordinate of the point in the center of a given small 
region of space. For x outside the area S the two states \Qx > and \Sx > are 
identical and therefore the overlap is unity. However for x inside S the states 
are different and the overlap is therefore some number e~'^ smaller than unity. 
The number of such regions inside the area is obviously of order S/P and we 
thus 

< UiS) >= exp{-7^} (54) 

In a weakly coupled theory this argument is confirmed by explicit calculation. 
The expectation value of the U-loop in the theory eq.(^) is given by the fol- 
lowing path integral 

< UiC) >= y^^dr exp|-y (a^0+z0x^)(a^0*-z0*XM)+A(0>-/i')'+c(0'^+(r) 

(55) 



with 



X^{x) = —Sf^'Sixo), xeS (56) 
= 0, x(^S (57) 



This expression directly follows from cq. (|5^) and integration over the canonical 
momentum in the phase space path integral. At weak coupling this path integral 
is dominated by a simple classical configuration. First, it is clear that the 
solution must be such that the phase of the field (p has a discontinuity of 2tt/N 
when crossing the surface S since otherwise the action is UV divergent due 
to singular x. Asymtotically at large distance from the surface the field should 
approach its vacuum expectation value. Since the source term x vanishes outside 
S, eveywhere where is continuous it has to solve classical equations of motion. 
Also, for values of xi and X2 which are well inside S the profile ((> should not 
depend on these coordinates, but should only depend on xq. It is easy to see 
that a solution with these properties exists: it is given by the "broken" domain 
wall solution. Recall that the vacuum is degenerate and so there certainly exists 
a classical solution of the equations of motion which interpolates between two 
adjacent vacuum states </) ^^o— >oo (j^o and </> —^xo—^-oo s^~4'o- Breaking this 
classical solution along the plane xq = and rotating the piece xq < by 27r/iV 
produces precisely the configuration with the correct boundary conditions and 
the discontinuity structure. The path integral in eq.(pq) is therefore dominated 
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by this classical configuration. Its action (up to corrections associated with the 
boundary effects of S) is aS where a is the classical wall tension of the domain 
wall which separates two adjacent Zjv vacua. Thus we find that the expectation 
value of the U-loop is related to the domain wall tension of the Zn domain wall 

by 

< U{S) >= exp{-aS} (58) 



4.2 U-loop in the unbroken phase 

Now consider the unbroken phase. Again the U-loop average has the form 



of the overlap of two states which factorizes as in cq.(53). Now however all 
observables noninvariant under vanish in the vacuum. The action of the 
symmetry generator does not affect the state |0 >. The state |S' > is therefore 
locally exactly the same as the state |0 > except along the boundary of the area 
S. Therefore the only regions of space which contribute to the overlap are those 
which lay within one correlation length from the boundary. Thus 

< UiS) >= cxp{-7P(5)} (59) 

where P{S) is the perimeter of the boundary of S. The absence of the area law 
is again easily verified by a perturbative calculation. In the unbroken phase the 
fiuctuations of the field (f) as well as the current density jo = iincj) — ■n*<j)*) are 
small. To leading order in the coupling constant 

< U{S) exp { - i / d^xSy < joix)My) > } (60) 

^ Jx,yeS ^ 

The possible area law contribution in the exponent is 



S d'x< joiO)joix) S lim G{p) (61) 
J p^o 

where G{p) is the Fourier transform of the charge density correlation function. 
The correlator of the charge densities however vanishes at zero momentum. 
This is because in the leading perturbative order the symmetry of the theory 
is actually U{1) and not just Zn as seen in eq. Since the vacuum state is 
invariant it follows that the total charge Q = J (fxjo{x) = jo{p = 0) on this 
state vanishes, and so does any correlation function that involves zero momen- 
tum component of the charge density. So the area contribution in eq.(pl) is 
zero. Strictly speaking in the leading order in perturbation theory eq.(|60) is 
not the complete result. The exact expression contains in the exponential also 
higher point correlators of the current density. Again however the possible area 
law contribution contains correlators of the total charge Q with powers of jo 
and therefore vanishes. 
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4.3 U-loop at high temperature 

Let us see now how the argument changes at high temperature. The important 
difference is that the vacuum is not a pure state but rather a statistical ensemble. 
The average of the U-loop is therefore not given by a single matrix element but 
rather by 

<U >^Y.e-^'^ <i\U\i> (62) 

i 

Let us consider the theory in which the Zn symmetry is broken at zero tem- 
perature. For concreteness we will think about Z2 symmetric theory, although 
qualitatively the discussion does not change for any TV. The two degenerate 
vacuum states are characterized by the value of the condensate < (j) >= i/^- 

In order to understand the behaviour of the U -loop we have to figure out what 
types of states contribute to the thermal ensemble. At zero temperature the 
only states that are of interst are those with finite energy. There are two towers 
of such states \n >^ and \n - constructed above each one of the degenerate 
vacua. These two towers of states do not talk to each other, not only because 
their overlap is zero, but also because they can not be connected to each other by 
action of any local (or semilocal) operator ^ < n|0|n' >-;i= 0. An immediate 
corollary of this is that a superposition of the type |a,/3 >= cy\n^ > 
violates clustering property of the correlators of local operators. For this reason 
at zero temperature in a spontaneously broken theory we are never interested 
in states which carry sharp quantum numbers of the broken symmetry. 

At finite temperature however we are also asking after states with finite energy 
density, and therefore infinite energy. This part of the spectrum looks rather 
different if the energy density involved is high enough. The two vacuum config- 
urations of the potential in eq.(|49|) are separated by a finite barrier. Let us call 
the hight of this barrier H. The states with energy density lower than H still 
separate into two towers. We will denote these states by \l >. However higher 
energy density states, with e > H have different nature. Their wave function is 
not localized in the field space to the vicinity of one of the vacua, but rather is 
spread over distances larger than the distance between the two vacuum states 
2/x. These states therefore naturally carry sharp quantum numbers with respect 
to the broken Z2 symmetry. These states we will denote by |/i >. In fact one 
expects that the higher the energy density the more these states look like the 
multiparticle states of a symmetric phase. That is to say as long as e >> it 
does not matter whether the potential has the double well structure or a single 
vacuum. These highly excited states should look like states with finite density 
of particles which carry the Z2 charge. 

At low temperatures, when the entropy effects are not important the contri- 
bution to the thermal ensemble comes only from the |^ > - sector since the 
Bolzman factor for any of the \h > states vanishes exponentially in the infinite 
volume limit. As we have argued earlier, the average of U in each one of these 
states has an area law behaviour and so does the whole temperature average of 
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u. 

When the temperature reaches Tc the phase transition occurs. The reason for 
the onset of the phase transition is that when the equihbrium energy density 
reaches critical threshold value, the \h > sector states start contributing to the 
thermal ensemble. The sudden change in the entropy due to these new channels 
drives the phase transition. Above the phase transition therefore there are two 
kinds of states that contribute to thermal averages. One can then write 

< U >T>T,= e^-^"/^ < n, l\U\n, I > + e^-^-^^ < n, h\U\n, h > (63) 

n s 

where n stands for all other quantum numbers. In fact once the entropy effects 
become important enough to excite the \h >-sector, the contribution of |Z > 
states to any physical observable becomes negligible. As discussed earlier each 
state in the second term gives a perimeter contribution exp{"fgP} to the average 
of the ?7-loop, so one could be tempted to conclude that the loop must have a 
perimeter law just like in zero temperature vacuum of a symmetric phase. This 
however is not necessarily the case. The reason is that < n, h\U\n, h > is not 
positive definite. In fact the number of states in which it is positive is roughly 
equal to the number of states on which it is negative. It is therefore very likely 
that the leading perimeter behaviour will cancel and the net result will again 
be an area law for < t/ >. 

Indeed if the ensemble can be thought of as an ensemble of Z2 charged free 
particles, the area law for [/-loop follows immediately]^ The [/-loop in such an 
ensemble is 

< U{S) >= J2 ^^^"(-1)" (64) 

Xi ,n 

where fi is the fugacity of a single particle, and the summation goes over all 
coordinates Xi of the particles inside the area 5* and over all possible numbers of 
particles n. Assuming that particles have a finite size A, so that there are S/A 
possibilities to place one particle inside the area S in the dilute gas approxima- 
tion the sum gives 

<UiS)>=exp{^^S} (65) 

We stress that the thermal ensemble of particles is Z2 invariant. The density 
matrix for such an ensemble can be written in the particle basis as 

n{x) 

The operator of the Z2 transformation acts on the n-particle states as 

U\n >= {-l)"\n > (67) 

^^This argument is borrowed from |20[.We thank Mike Teper and Biaggio Lucini for discus- 
sions of this point. 
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and so 



UpU^ = p 



(68) 



The explicit simple formula eq.(|65|) is derived in the dilute gas approximation. 
We expect that the physics will be similar as long as the interaction between the 
particles is short range. Whenever the interaction is long range the behaviour of 
< U{S) > can be different. For instance one does not expect area law behaviour 
if the particles are bound into pairs since in this case only Z2 invariant states 
contribute to the thermal ensemble. 

Our conclusion is that at finite temperature the behaviour of the ?7-loop is not 
strongly related to the mode of the realization of the symmetry. It is rather 
more likely to have an area behaviour. 

To reiterate, the physics involved is very simple. At zero temperature when 
acting on a state, the [/-loop performs the Zn transformation inside the loop. 
The only degrees of freedom that are changed by this operation inside the loop, 
are the Z^ - noninvariant fields. If the vacuum wavefunction depends on the 
configuration of the noninvariant degrees of freedom (the state in question is not 
Z^ invariant) the action of [/-loop affects the state everywhere inside the loop. 
The VEV of [/-loop then falls off as an area. If the vacuum is Zjv invariant, the 
wavefunction does not depend on the configuration of the noninvariant degrees of 
freedom. The action of U -loop then perturbs the state only along the perimeter, 
hence the perimeter law in the unbroken phase. 

At finite temperature however the thermal ensemble may even in the symmet- 
ric phase contain significant contributions from states with nonvanishing Z^ 
charges. The [/-loop therefore perturbs the thermal ensemble very significantly 
everywhere inside the area, and the natural outcome is an area law. 

The argument is quite general and does not depend on the exact form of the 
Zn invariant potential and more generally on the field content of the theory - 
we could have added any number of extra fields to the theory eq.(|49|) without 
changing the conclusions. The same relation must exist between the mode of 
realization of the magnetic Z^ symmetry and the behaviour of the Wilson loop 
in the pure Yang-Mills theory. The direct analogs of the scalar field 0(a;) in 
eq.(^) and the U-loop of the scalar theory are correspondingly the vortex field 
V{x), and the spatial Wilson loop W{C). 

As we have shown in the previous section, the magnetic Zn is restored at high 
temperature. The Wilson loop is nevertheless likely to have an area law as 
is indeed indicated by all existing lattice data. In this context we note that 
analytic strong coupling results also give area behaviour [p|. 
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4.4 Wilson loop in 3+1 dimensions 



The previous considerations generalize to the 3+1 dimensons. At zero tem- 
perature in the broken phase when acting with the Wilson loop W{C) on the 
vacuum one changes the state of those magnetic vortices which loop through C. 
The number of such vortices which are present in a generic configuration in the 
broken phase is proportional to the minimal area subtending C. The number of 
the degrees of freedom that is changed by the action of W is thus proportional 
to the area S. Each of these degrees of freedom contributes a factor smaller 
than unity to the overlap with the vacuum state and so the VEV of W scales 
with the exponential of the area. In the unbroken phase the vacuum does not 
contain vortices of arbitrarily large size. The size of the vortices present in the 
vacuum is cutoff by the relevant correlation length. This is the case if the gauge 
group is completely broken by the Higgs mechanism. Therefore for contours 
C of linear dimension much larger than this length, the action of W{C) only 
disturbs degrees of freedom close to the contour C itself and the VEV must have 
the perimeter behaviour. 

At high temperature even the symmetric thermal ensemble is populated by 
vortices. These vortices are not free since apart from them the ensemble also 
contains "free" charges. However unless this background of charges induces long 
range interactions between the vortices, the most probable result for the Wilson 
loop is the area law. A more detailed knowledge of vortex dynamics is necessary 
to draw a firm conclusion. 

To close this section we note that the present considerations do not apply to 
Abelian theories. The magnetic symmetry does exist in this case too, but here it 
is the continuos U (1) group and the spectrum is massless. In this case there is no 
reason to expect the local factorization of the overlap and generically therefore 
the arguments of this section do not hold. In particular in the presence of long 
range correlations it is perfectly possible that the Wilson loop has a perimeter 
law even though the state is perturbed everywhere inside the area boundede by 
the loop[3. 

5 Discussion. 

In this paper our aim was to point out two facts. First that the calculation of 
the VEV of the 't Hooft loop 0] implies the restoration of the magnetic Zn 
symmetry above the deconfinement transition. Second, that the mode of real- 
ization of magnetic symmetry is closely related to the behaviour of the spatial 
Wilson loop. 

^^In 2+1 dimensions it is actualy only the noncompact Abelian theories that are excluded 
from the consideration. Compact theories are massive and therefore should behave in the 
same way as the nonabelian Yang-Mills. 



26 



At zero temperature this relation is very rigid: spontaneous breaking of Zj^ 
implies the area law behaviour for W while unbroken Z^r leads to perimeter 
law behaviour. At high temperature however even though the Zj^ symmetry is 
restored, the Wilson loop may have an area law. This is the consequence of the 
fact that even a Zjq - invariant the thermal ensemble can contain a significant 
contribution of Z]^ nonsinglet states. The area law is particularly simply un- 
derstood if the thermal ensemble at high temperature is well appproximated by 
an ensemble of weakly interacting magnetic vortices. The vortex gas argument 
has been previously brought up in favour of the area behaviour of the spatial 
Wilson loop in This behaviour is also confirmed by several lattice gauge 

theory calculations |l6|| . 

An alternative possibility is that due to the as yet unknown vortex dynamics, 
there is vortex - antivortex binding. To explore such a possibility it would be 
very interesting to measure on the lattice the free energy of a magnetic vortex. 
In ref.(|]l^) the behaviour of the free energy of magnetic and electric fluxes has 
been discussed in the low temperature phase. To be able to do it in the lattice 
framework one has to define the theory in a finite volume. As discussed by 't 
Hooft this can be achieved by imposing on the potentials periodic boundary 
conditions modulo a gauge transformation. As discussed in ref. jist this admits 
the presence of vortices in 2+1 and of the vortex lines in 3+1 dimensions, 't 
Hooft's discussion was based on a Euclidean rotation identity for the twisted 
4d path integrals valid for any temperature, and a factorization property of 
magnetic and electric fluxes. In the notation of ref.[|l5j: 

F{e,m)^F,{e)+F^{m) (69) 

Its validity at low T is very reasonable, but is inconsistent with the Euclidean 
rotation identity at high T. Based on this 't Hooft could prove (N < 3) that 
in the confining phase, where the free energy of an electric flux is linear with 
the length (with the string tension p), the free energy of magnetic flux vanishes 
exponentially in the infinite volume limit. For a magnetic flux in, say the z- 
direction it is exp ^pLxLy. Thus the free energy of a magnetic flux is related 
to the behaviour of the Wilson loop. The free energy of an electric flux in the z 
direction in the hot phase vanishes exponentially like exp —aLxLy where a is the 
surface tension found in ref. |^ . So the next obvious question is how the magnetic 
flux free energy behaves in the hot phase. In a vortex gas picture this free energy 
should vanish exponentially in the infinite volume limit. Such a calculation 
in 2+1 dimensions has been performed and (modulo some uncertainty related 
to imperfect measurement of global vorticity) results are consistent with this 
expectation It would also be instructive to see how in the hot phase the 
additivity of electric and magnetic fluxes is broken. 

We note that a recent lattice calculation |2^ measures the monopole-antimonopole 
correlation. The results of ||2l| point to the screened behaviour of this corre- 
lation function for all temperatures. So in the hot phase it behaves like its 
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electric partner, the correlator of Polyakov loopq^. This via 't Hoofts argu- 
ment, is consistent with the measured area behaviour of Wilson loops and 
would imply that the magnetic flux free energy would fall off with an area law 
for all temperatures. 

It is interesting to note that the vortex gas picture is equally applicable in 
high temperature confining and nonconfining gauge theory. In particular one 
can consider an SU (N) gauge theory with sufficient number of adjoint Higgs 
fields, so that the gauge theory at zero temperature is broken completely. In 
this situation the magnetic Zn symmetry is unbroken in the vacuum and the 
Wilson loop has a perimeter law. Magnetic vortices are finite energy excita- 
tions with the mass of order M = M'^/g^, where is the vector boson mass. 
When the system is heated one expects that the thermal ensemble will con- 
tain a dilute gas of these vortices at any temperature. Therefore at any finite 
temperature the spatial string tension should be nonzero, although at low tem- 
peratures it will be exponentially suppressed if the theory is weakly coupled: 
aocM2exp{-M2/52r}. 
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